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Abstract 

Non-Hermitian but 'P^'T^-symmetrized spherically-separable Dirac and 
Schrodinger Hamiltonians are considered. It is observed that the descen- 
dant Hamiltonians H r , He, and H v play essential roles and offer some 
"user-feriendly" options as to which one (or ones) of them is (or are) non- 
Hermitian. Considering a P^T^-symmetrized H v , we have shown that 
the conventional Dirac (relativistic) and Schrodinger (non-relativistic) en- 
ergy eigenvalues are recoverable. We have also witnessed an unavoidable 
change in the azimuthal part of the general wavefunction. Moreover, set- 
ting a possible interaction V (8) 7^ in the descendant Hamiltonian He 
would manifest a change in the angular ^-dependent part of the general 
solution too. Whilst some "P^T^-symmetrized H v Hamiltonians are con- 
sidered, a recipe to keep the regular magnetic quantum number m, as 
defined in the regular traditional Hermitian settings, is suggested. Hamil- 
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tomans possess properties similar to the "PT-symmetric ones (here the 
non-Hermitian P^T^-symmetric Hamiltonians) are nicknamed as pseudo- 
VT -symmetric. 

PACS codes: 03.65.Ge, 03.65.Ca 

Keywords: Non-Hermitian Hamiltonians, spherical-separablility, pseudo- 
PT-symmetry. 

1 Introduction 

In the search for the reality conditions on the energy spectra/eigenvalues of 
non-Hermitian Hamiltonians [1-35], it is nowadays advocated (with no doubts) 
that the orthodoxal mathematical Hermiticity requirement to ensure the re- 
ality of the spectrum of a Hamiltonian is not only fragile but also physically 
deemed remote, obscure and strongly unnecessary. A tentative weakening of 
the Hermiticity condition through Bender's and Boettcher's [1] PT-symmetric 
quantum mechanics (PTQM) (with V denoting parity and T time-reversal, a 
Hamiltonian H is PT-symmetric if it satisfies VTHVT = H) has offered an 
alternative axiom that allows for the possibility of non-Hermitian Hamiltonians. 

Such a PTQM theory, nevertheless, has inspired intensive research on the 
non-Hermitian Hamiltonians and led to the so-called pseudo-Hermitian Hamil- 
tonians (i.e., a pseudo-Hermitian Hamiltonian H satisfies rj H rj^ 1 = or 
r\ H = r/, where rj is a Hermitian invertiblc linear operator and (t) denotes the 
adjoint) by Mostafazadeh [16-21] which form a broader class of non-Hermitian 
Hamiltonians with real spectra and encloses within those PT-symmetric ones. 
Moreover, not restricting r\ to be Hermitian (cf., e.g., Bagchi and Quesne [33]), 
and linear and/or invertible (cf., e.g., Solombrino [28], Fityo [29], and Mustafa 
and Mazharimousavi [30-32]) would weaken pseudo-Hermiticity and lead to real 
spectra. 
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Based on the inspiring example, nevertheless, by Bender, Brody and Jones 
[35] that H = p 2 + x 2 + 2x is a non-'PT-symmetric whereas a simple amendment 
H = p 2 + (x + l) 2 — 1 (that leaves the Hamiltonian invariant and allows parity to 
perform reflection about x = — 1 rather than x — 0) would consequently classify 
H as PT-symmetric (i.e., reflection need not necessarily be through the origin) 
and promoting ZnojiPs understanding [34] of Bender's and Boettchcr's PTQM 
(i.e., V and T need not necessarily mean just the parity and time reversal, 
respectively), we may introduce [36] a time-reversal- like, 

T v : L z — -id/dip — ► -L z = id/dip, ip — ► ip,i — ► —i (1) 

and a parity-like 

: L z — ► -L„ <p — » (2tt - ip) , (2) 

operators that might very well be accommodated by Bender's and Boettcher's 
PTQM. 

In this case, V v acting on a function / (r, 6, ip) e L 2 would read P v f (r, 6, ip) = 
/ (r, 6, 2n — ip). Moreover, / (r, 9, <p) is said to be ^7^,-symmetric if it satis- 
fies V^Tpf (r,6,<p) = f(r,6,<p). Hence, our new operators leave the coordi- 
nates r and 9 unaffected and are designed to operate only on the azimuthal 
descendent eigenvalue equation (e.g., equation (21) below) of the spherically- 
separable non-Hermitian Hamiltonians (4) and (5). However, it should be 
noted that our parity-like operator P v in (2) is Hermitian, unitary, and per- 
forms reflection through a 2D-mirror represented by the xz-plane. Yet, the 
proof of the reality of the eigenvalues of a 'P ¥ ,7^-symmetric Hamiltonian is 
straightforward. Let the eigenvalue equation of our 'P v 7^-symmetric Hamil- 
tonian be Hip (r, 6, <p) = Eip (r, 6, ip), then P^T^Hip = P v T v Eip = Eip. Using 
[P V T V ,H] =0we obtain Eip — E*ip and E is therefore pure real (in analogy 
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with Bender, Brody and Jones in [35] and fits into PTQM-recipe). 

In the forthcoming proposal, using spherical coordinates, we depart from 
the traditional radial potential setting (i.e., V (r) = V (r)) into a more general 
potential of the the form 



V (r) = V (r, 0,ip)=V (r) + 



V{0) + V(<p) 



(3) 



r 2 sin 2 9 

We shall use such potential setting in the context of Schrodinger Hamiltonian 

H = -V 2 +V(r), (4) 

and within an equally-mixed vector, V (r), and scalar, S (r), potentials' setting 
in the Dirac Hamiltonian 

H = a-p + p[M + S(r)] + V{r) 1 (5) 

with the possibility of non-Hcrmitian interactions' settings in the process. How- 
ever, it should be noted that such interactions in (3) with V (r) = —a/r, 
V (9) = —b 2 , and V ((f) = represent just variants of the well known Hart- 
mann potential [38-47] used in the studies of ring-shaped organic molecules. 

For the sake of making our current proposal self-contained, we revisit, in 
section 2, Dirac equation in spherical coordinates and give preliminary foun- 
dation on its separability. We connect, in the same section, Dirac descendant 
Hamiltonians with those of Schrodinger and provide a clear map for that. In 
section 3, we explore some consequences of a class of complexified but V V T V - 
symmetrized azimuthal Hamiltonians. For a complexified azimuthal interac- 
tion V (ip) € C (with V (r) , V (6) € R) we use three illustrative examples for 
V(9) = 0, V{6) = 1/2, and V (9) = l/(2cos 2 6>). In section 4, a recipe of 
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generating functions is provided to keep the magnetic quantum number as is, 
whenever deemed necessary of course. In the process of preserving the magnetic 
quantum number to, a set of isospectral 92-dependent potentials, V ((f), for each 
set of V (r) and V (9) is obtained. This would, moreover, allow reproduction of 
the conventional-Hcrmitian rclativistic and non-relativistic quantum mechani- 
cal eigenvalues within our 7^7^-symmetric non-Hcrmitian settings. We give 
our concluding remarks in section 5. 

2 Separability and preliminaries of Dirac and 
Schrodinger equations revisited 

Dirac equation with scalar and vector potentials, S (r) and V(r), respectively, 
reads (in H = c = 1 units) 



{a • p + (3 [M + S (r)] + V (r)} V (r) - (r) , 



(6) 



where 



p = — iV , a = 



er 



<T 



I 
-I 



(7) 



and a is the vector Pauli spin matrix. A Pauli-Dirac representation would, with 

/ 

V 



V(r) = 




(8) 



yield the decoupled equations 



(a ■ p) x 2 (r) = [E — V (r) — M — S (r)] Xl (r) , 



(9) 



(a ■ p) xi (r) = [E - V (r) + M + S (r)] X2 (r) . 



(10) 
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An equally-mixed scalar and vector potentials (i.e., S (r) = V (r)) leads to 



X2 (r) = 



g ■ P 

E + M' 



(11) 



and 



[- V 2 + 2 (£ + M) V (r)] Xi (r) - [£ 2 - M 2 ] Xl (r) . 



(12) 



Departing from the traditional "just-radially-symmctric" vector potential (i.e., 
V (r) = V (r)) into a more general, though rather informative, vector potential 
(in the 3D spherical coordinates r, 0, and <p) of the form 



V (r) = V (r, 0, = V (r) + 



r 2 sin 2 I 



(13) 



would, with 



Xi (r)=X! (r,fl >¥ ))=iJ(r)e(fl)*( V ) ; 



(14) 



imply 



R(r) 



1 {^{ r2 i)- 2{E+ M) v {r) r2 + {e2 m2 ) r2 } R (r) 

1 r d ( . „ 9 \ 2(E + M)V(e)~\ _ ... 

B (0) 



6(0) sin L<90 V 30 



d ( . n d\ 2(E + M) V (9) 
sin ! 



sin0 



+ 



$ (y>) sin 2 \dip 



d 2 

— - 2 (E + M) V {ip) ) $ (p) = 



The separability of which is obvious and mandates 



{ r 2 dr ( 



2 -^j- ^ - 2 (E + M) V (r) + (£ 2 - Af 2 ) J R (r) = 0, 



1 d ( . d_ 
sm9 d9 \ Sm °d6 



2 (E + M) V (0) + m? 



sin 



A 



6 (9) = 0, 



- 2 (£ + M ) V (ip) + m 2 ) $ (p) = 



(15) 



(16) 

(17) 
(18) 
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where m 2 and A are separation constants to be determined below. Yet, in a 
straightforward manner, one can show that both Dirac and Schrodinger equa- 
tions (with V (r) in (13)) would read 



where 



1 d ( 2 d \ A 



^fr( r Tr)-^- V °"^ +X } R ^=°> 



1 d ( . d_ 
sine dO \ d6 



d 



V eff (9)+m 2 
sin 2 9 



A 



e (0) = o, 



^2-K//(v)+m^)*(^)=0 



Veff (r) = < 



Veff (0) 



Veff {f) = < 



A = <^ 



V (r) for Schrodinger 
2 (E + M) V (r) for Dirac 

V (9) for Schrodinger 
2 (E + M) V (9) for Dirac 

V (ip) for Schrodinger 
2 (E + M) V (ip) for Dirac 

E for Schrodinger 

E 2 - M 2 for Dirac 



(19) 

(20) 
(21) 

(22) 

(23) 

(24) 

(25) 



The map between Schrodinger and Dirac equations is clear, therefore. Moreover, 
one can safely name three "new" descendant Hamiltonians and recast the cor- 
responding eigenvalue equations (with A = E for Schrodinger and A = E 2 — M 2 
for Dirac) as 

H r R (r) = XR (r) , H e Q (9) - A9 (9) , {ip) = m 2 $ (<p) . (26) 



Of course it is a straightforward to work out the explicit forms of H r , Hq, 
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and Hp from (19), (20), and (21), respectively. Moreover, if we substitute 
U(r) = R (r) /r in (19) then U (0) = = U (oo). Yet, whilst 6 (0) and 6 (tt) 
should be finite, $ (<p) should satisfy the single-valuedness condition $ (ip) = 
$(</? + 27r). At this point, we argue that the reality of the spectrum of Dirac 
eigenvalue equation (6) is ensured not only by requiring m, A, A e R but also by 
requiring K 9 A + M 2 = E 2 > 0. With this understanding, we may now seek 
some 'PT-symmctrization recipe (be it Lcvai's [37] regular 'PT-symmetrization 
or "P^T^-symmetrization of Mazharimousavi [36]) for each (at a time) of the 
descendant Hamiltonians in (26). 

3 Consequences of complexified 7^7^-symmetrized 
azimuthal Hamiltonians 

The eigenvalue equation in (21) with a € M as a coupling parameter in a 
complexified-azimuthal effective interaction of the form 

Veff M - -a 2 e* (27) 



would read 

d 2 



d(p z 



$ (ip) = 0. (28) 
Hence, a change of variable of the form z = e %ip l 2 would result in 

a ^*W +z «(£)_ (4ma + 4fl V)* W = 0. (29) 



z 



dz 2 dz 



Obviously, equation (29) is the modified Bessel equation with imaginary argu- 
ment and has two independent solutions. The linear combination of which reads 
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the general solution 

$ (z) = dhm (2az) + C 2 K 2m (2az) . 

Each of these independent solutions should identically satisfy the single- valuedness 
condition $ (<p) = <& (<p + 2ir) . One may, nevertheless, use the identities of [48] 
and closely follow Mazharimousavi's treatment (namely, equations (17) - (28) 
in [36]) and show that 

I 2m (2ae^+ 2 ^/ 2 ) = I 2m (2ae^' 2 ) , (30) 

K 2m (2ae J ^+ 27r )/ 2 ) ^ K 2m (2ae^ 2 ) . (31) 
Therefore, the regular solution collapses into 

$ (z) = C x I 2m (2az) = C x I 2m (2ae^/ 2 ) , (32) 

Under such settings, it is obvious that the Hamiltonian represented in (28) 
reads 

H v = -^-W*, (33) 
and qualifies to be a T-^T^-symmetric non-Hermitian Hamiltonian. That is, 

V T H V T =77 

' tp •*• if J J if ' <p If — ±J -ip- 

On the other hand, the eigenvalue equation (21) with H v would admit a regular 
azimuthal solution represented by the modified Bcsscl function 

* (</>) = C m . a I 2m (2ae^/ 2 ) ; m = 0, ±1, ±2, ■ ■ ■ . (34) 
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and satisfies the single-valuedness condition $ (ip) = $(</? + 2n) with C m . a as 
the normalization constant to be found through the relation 



Hereby, we have used the fact that our $ m (ip) in (34) is T^^-symmetric sat- 
isfying V v T v <£ m (ip) = $ m (ip). 

This would, in effect, suggest that since R 3 m = 0, ±1,±2, ■ ■ ■ and Hg of 
(20) is therefore Hcrmitian, then Hg of (20) admits real eigenvalues represented 
by A e R. Some illustrative consequences (with Hg of (20) kept Hermitian) are 
in order. 

3.1 Consequences of V (9) = in (23) 

Should V (9) = 0, one may clearly observe that equation (20) is the very well 
known associated Legendre equation in which A = 1(1 + 1), where I is the 
angular momentum quantum number, and 9 (8) = P™(cos#) are the associated 
Legendre functions. Hence, following the regular textbook procedure one may, 
in a straightforward manner, show that I > |m| (i.e., m = 0, ±1, ±2, • • • , ±£, is 
the regular magnetic quantum number). 

Consequently, as long as the Hermitian radial equation (19) is solvable (could 
be exactly-, quasi-exactly-, conditionally-exactly-solvable, etc) for the radial 
interaction V e ff (r), the spectrum remains invariant and real. However, the 
global wavefunction 



(with n r = 0, 1, 2, • • • as the radial quantum number) would indulge some new 



l = ($ m (ip) /V v T v $ m (p)) = |C, 



/ " hrahae^' 2 ) 2 dip. (35) 



Jo 



Xi (r,0,<p) = ip S ch {r,0,<p) 
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probabilistic interpretations. This is due to the replacement of the regular 
spherical harmonics Yi m (0, ip) part (for the radially symmetric 3D-Hamitonians) 
by the new T-^T^-symmetric part P™(cos 9)& m (p) (defined above for our V V T V - 
symmetric non-Hcrmitian Hamiltonian model). 

To see the effect of such a "P v 7^-symmetrization on the probability density, 
we consider a radial Coulombic effective interaction V e ff ( r ) = — 1/f accom- 
panied by an azimuthal effective interaction V e ff (p) = — Cretan illustrative 
example of fundamental nature). In figures 1 and 2 we plot the correspond- 
ing probability densities at different values of the coupling parameter a for the 
principle quantum numbers n = 1 and n = 2 for I = = m. It is clearly 
observed that whilst the probability density for small a imitates the Hermi- 
tian 93-independent probability density trends, it shifts and intensifies about 
\<p>\ = as a increases (indicating that the corresponding state is more localized, 
therefore). In this case, of course, the rotational symmetry of a purely "just- 
radially-symmetric" Coulombic interaction breaks down as a result of V e fj (ip) . 



3.2 Consequences of V (9) = 1/2 or V (6) = 1/(2 cos 2 #) in 



where rh = \[E + M + m 2 for Dirac and fh = y/l/2 + m? for Schrodinger 
settings. Similar equation was reported by Dutra and Hott [47]. The regular 
solution of which can (taking a — f3 — and 7 = 1 in equation (12) of ref. [47] 
to match with our settings) very well be copied and pasted to read (for Dirac 



(23) 



Taking V (6) = 1/2 in (23) would imply 




(37) 
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equation) 



0(6) = yr(l-y) v 2 F 1 (-k,b;d;y) ; y = cos 2 (9/2) . 



(38) 



where 



v = p 



\\fm 2 + E + M ; b = k + 4v + 1; d = 1 + 2v, 



fc = 0, 1, 2, ■ ■ ■ is a "new" quantum number, and 



(39) 



A = \ (b + kf 1 = 1 [(b + fc + 1) (6 + fc - 1)] . 



(40) 



On the other hand, V (0) = 1/ (2 cos 2 0) would (taking a = /3 = 0and7 = l 
in equation (13) of ref.[12) for Dirac equation) result in 



P=\ + ^l + A{E + M)- v = ^y/m?+E + M 



& = fc + 2(p + w) + i; d=2p+i 



and 



A= (b + kf 



1 



b + k+l) (b + k-\ 



(41) 
(42) 

(43) 



For Schrodinger case, nevertheless, one may just replace the term (E + M) by 
(1/2) in the above expressions and get the corresponding eigenvalue results. 
Then the general solution for both cases would read 

Xi {r, 0, <p) = i>sch (r, 8, V>) 

= N nT ^ m R nT . k (r) yP (1 - y) v 2 F 1 (-fc, 6; d; y) I 2m (2ae^ 2 ) , (44) 



where N nri k, m is the normalization constant that can be obtained in a straight- 
forward textbook procedure. Hereby, we witness that the general solution (44) 
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exhibits the change not only in the azimuthal part but also in the angular 6*-part. 



4 Preservation of the magnetic quantum num- 
ber m and isospectrality 

To keep the magnetic quantum number as is (i.e., m — 0, ±1, ±2, • • • ), one may 
consider the azimuthal part of the general solution to be of the form 

<f m (ip) = e im <?F (ip) , (45) 

where F (ip) satisfies the single- valuedness condition F (ip) = F (ip + 2n). 

Under such setting, the corresponding eigenvalue equation in (21) (with 
primes denoting derivatives with respect to ip) reads 

F" (<p) + 2imF' (ip) - V ef f (ip) F (p) = 0. (46) 

In this case F (ip) would serve as a generating function for the sought after 
azimuthal potential V e fj (p) and shapes the form of the azimuthal solution 
$ m (ip). As an illustrative example, a generating function F (p) — cos<p would 
imply 

v ef f (v) = - [1 + 2im tan <fi] ( 47 ) 

which is indeed a non-Hermitian and T^T^-symmetric, V v T v V e ff (<p) — V e f / (p) . 

However, one may wish to follow the other way around and consider a V V T V - 
symmetric V e /f (<p) and solve (46) for F (p). In this manner, V e j f ('p) would now 
serve as a generating function for F (p) and consequently a generating function 
for $ m (ip). An immediate example is in order. Consider 

Veff(f) = ~e^ (48) 
13 



and solve (44) for a regular F (ip) to obtain 



R3m = 0,il,i2,-. 



(49) 



Then, (45) would read 




(50) 



It is, therefore, obvious that all effective potentials V e // (<p) satisfying (48) would 
essentially change the azimuthal part of the general solution. 

Moreover, in the process of preserving the magnetic quantum number m 
as defined in the regular Hermitian settings, a set of isospectral (^-dependent 
potentials, V e ff (tp), is obtained. That is, for each set of V e /f (r) , V e /f (9) £ K, 
all (^-dependent potentials, V e fj (ip), satisfying (46) are isospectral. 

5 Concluding remarks 

In the build up of a generalized quantum recipe (Bender's and Bocttcher's 
PTQM in this case), a question of delicate nature arises in the process as to 
"would PTQM be able to recover some results (if not all, to be classified as 
a promising theory) of the conventional Hermitian quantum mechanics?". To 
the best of our knowledge, only rarely and mainly within regular Hermitian 
(but PT-symmetric) settings examples were provided such as the one by Ben- 
der, Brody and Jones [35] mentioned in our introduction section above (i.e., 
H = p 2 + x 2 + 2x) . The reality of the energy eigenvalues and other quantum 
mechanical properties (rather than the " recoverability of Hermitian quantum 
mechanical" results) were the main constituents and focal points in the studies 
of the non-Hermitian 'PT-symmetric Hamiltonians. In our current proposal, 
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with a new class of non-Hermitian 7^7^-symmetric Hamiltonians (having real 
spectra identical to their Hermitian partner Hamiltonians), we tried to fill this 
gap, at least partially. 

Through our over simplified non-Hermitian T-^T^-symmetrized Hamiltonian 
(33), we have shown that some conventional relativistic and non-relativistic 
quantum mechanical results are indeed recoverable (the energy eigenvalues here). 
We have witnessed, however, an unavoidable change in the azimuthal part of 
the general wavefunction. Such a change would introduce some new probabilis- 
tic interpretations. With V (0) =0 and V e ff (<p) = —a 2 e lip , for example, we 
have observed that a quantum state becomes more localized as the probability 
density intensifies at a specific point \ip\ = (documented in figures 1 and 2 
of section 3.1). Moreover, setting V (0) ^ (again with V e ff (ip) — —a 2 e tip in 
H v ) in the descendant Hamiltonian Hg has indeed manifested a change in the 
angular ^-dependent part of the general solution too (documented in section 
3.2). This would, of course, has some "new" effects on the probabilistic inter- 
pretations in turn. Yet, a recipe to keep the magnetic quantum number m as 
defined in the regular Hermitian quantum mechanical settings is suggested. 

In connection with the current proposal's spherical-separability and non- 
Hcrmiticity, it is obvious that the descendant Hamiltonians H r , Hg, and H v 
play essential roles and offer some " user- friendly" , say, options as to which 
one (or ones) of them is (or arc) non-Hermitian. Be it ^7^-symmetric, PT- 
symmetric, pseudo-Hermitian or 77-pseudo-Hcrmitian, they very well fit into 
Bender's and Boettcher's PTQM (irrespective with their nicknames and with 
the understanding that V and T need not necessarily identify just parity and 
time reversal, respectively). Yet, a complexification of 7^ V (0) € C in Hg 
with the understanding that a parity-like Vg and a time reversal-like Tg oper- 
ators may very well suggest a similar "PeTg-symmetric Hg Hamiltonian. Such 
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non-Hermitian ^7^-symmetrized and/or Te% -symmetrized (anticipated to be 
feasible but yet to be identified) Hamiltonians better be nicknamed as pseudo- 
VT -symmetric Hamiltonians. 
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Figures' Captions: 

Figure 1: Shows the effect of V e ff (<£>) ; 
as the coupling parameter a increases for n 

Figure 2: Shows the effect of V e ff (f) z 
as the coupling parameter a increases for n 



-a 2 e l(p on the probability density 

1, £ = = m. 

-a 2 e l(p on the probability density 

2, £ = = m. 
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